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ABSTRACT 

We examine the ability of the Sunyaev-Zel'dovich effect (SZE) statistics pro- 
posed by Lee (2002) using numerical simulations. The statistics describe the 
distribution of the peak heights in the noise-contaminated SZE sky maps, and 
provide an analysis technique to sort out the noise-contamination effect and es- 
timate the number density of real clusters. The method is devised to be suit- 
able for the interferometric SZE observations in drift- scanning mode like AMiBA 
experiment. We apply the proposed method to a set of realistic SZE sky maps 
constructed from large-scale cosmological simulations, and show that the method 
indeed allows us to estimate the number density of clusters efficiently. The ef- 
ficiency of the method is demonstrated in two aspects: (i) it can reconstruct 
the number density of clusters even though the clusters contribute only a small 
fraction (~ 10%) of the total number of peaks in each SZE map; (ii) it can count 
even those clusters whose amplitudes are low enough to be compatible with the 
noise level while conventional method of cluster identification using the threshold 
cutoff cannot count them properly within the same observation time. Thus, the 
proposed method may be useful for the study of cluster abundance exploiting 
future interferometric SZE surveys such as AMiBA. 

Subject headings: cosmology: theory — SZ effect: clusters 



1. INTRODUCTION 

The thermal Sunyaev-Zel'dovich effect (SZE) by galaxy clusters represents the spec- 
tral distortion of the cosmic microwave background radiation (CMBR) through the inverse 
Compton scattering by the hot intra-cluster electrons (Sunyaev & Zel'dovich 1972). Since 
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the SZE is independent of the redshift, it provides a unique tool to detect high-redshift clus- 
ters. Thus, the blank-field SZE survey can measure the formation epoch and the abundance 
evolution of galaxy clusters, which can be used as a probe to put constraints on the cosmo- 
logical parameters (Barbosa et al. 1996; Henry 1997; Bahcall & Fan 1998; Viana & Liddle 
1999; Borgani et al. 1999; Fan & Chiueh 2001; Grego et al. 2001; Molnar et al. 2002). 

The SZE flux diminution is directly proportional to the cluster Comptonization param- 
eter y, defined as 

y0) f ne{W)T,{W)dl, (1) 

where 9 represents the location vector along the line-of-sight {9 = 9/\9\) and the integral is 
over the CMBR path in the line-of-sight direction. Hence the SZE survey at a fixed frequency 
basically provides a map of the ^/-parameter. The number density of galaxy clusters, however, 
cannot be directly measured from the observed SZE maps because of the presence of noise. 

Conventionally, signals from the noise-contaminated maps are identified as the local 
peaks whose flux amplitudes (for the SZE observation, proportional to the peak height of the 
Comptonization parameter, y) signiflcantly exceed the noise level. Therefore, to detect lower- 
amplitude signals using the conventional approach, one has to decrease the noise level by 
increasing the observational time. Unfortunately, the necessary observational time increases 
inversely proportional to the squares of the imposed noise level. It is clearly inefficient to 
use the conventional signal-selection method for the estimate of the cluster number counts, 
because of the high costs of the interferometric observations. 

Recently Lee (2002) proposed a statistical strategy to estimate efficiently the number 
density of galaxy clusters from the noise-contaminated SZE maps constructed from the in- 
terferometric observation in drift-scanning mode, with the AMiBA experiment as a target 
model. By the efficient estimate, we mean that the method is capable of estimating the 
number density of clusters which contribute only a small fraction of the total number den- 
sity of peaks in the SZE maps, counting properly even those clusters whose peak heights are 
low enough to be compatible with the noise level without increasing the observational time. 
The key feature of the method is to exploit the property that the noise-contamination is de- 
scribed by a Gaussian statistic. Then the cluster number density can be estimated by simply 
deconvolving a one-dimensional Gaussian function and the zero-th order approximation to 
the cluster number density. Although the method has been shown to render reasonable re- 
sults for idealized Monte-Carlo simulations (Lee 2002), it is necessary to test its ability by 
applying it to more realistic simulations before using it in real SZE surveys. 

In this Letter, we test the ability of the SZE statistical strategy using the outputs of a 
large cosmological simulation. In §2, we provide a concise overview of the model experiment 
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AMiBA and the SZE statistics. In §3, we describe the numerical simulation, explain how- 
to apply the statistical strategy to the simulated SZE maps, and compare the reconstructed 
cluster number densities with the original ones. In §4, we summarize the results, and give 
concluding remarks. 



2. The SZE CLUSTER SURVEY 

2.1. AMiBA in a Drift-Scan Mode 

For the blank-field SZE survey, we have two types of observations: the single dish 
bolometric observations such as BOLOCAM (Glenn et al. 1998) and the multiple dish in- 
terferometric observations such as our model experiment, AMiBA (Lo et al. 2000). AMiBA 
(Array for Microwave Background Anisotropy) is an interferometer telescope with 19 dishes 
of 1.2 meter, dedicated to SZE observation at operating frequency of 90 ± 16 GHz. It is 
designed to be most suitable for a deep survey of low-mass clusters at high redshifts. For 
this purpose, it has relatively small beam of 2 arcminute, and plans to survey 5 square degree 
for 20 hours per field on average (the total period is over 7 months), expected to be able to 
detect 1 cluster per every 9 hour. 

An observational strategy of AMiBA is to use the drift-scan method, which will help 
keep the ground fringes constant with time, improve the Fourier-space resolution, and most 
importantly simplify the noise analysis (Pen et al. 2002). A drift-scanned sky map from 
AMiBA experiment will be a combination of a small fraction of cluster signals (~ 10%) and 
a large fraction of dominant instrumental noise. The distribution of cluster signals is non- 
Gaussian, while that of the instrumental noise is Gaussian with no spatial correlation. Since 
the statistical property of a Gaussian field is well known, the noise analysis for the SZE sky 
maps from AMiBA experiment in drift-scanning mode may be analytically tractable (Zhang 
et al. 2002; Lee 2002). 



2.2. The SZE Statistics : An Overview 

The peak number density of a Gaussian random field is analytically given by (Longuet- 
Higgins 1957; Bond & Efstathiou 1987): 

( \ 1 0-2 -4 Tr 2^ exp[-l(x-7//)V(l-7^)] ^ ... 
(^^* ^ " ^ [27r(l-7^)]V2 (2) 
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where v is the peak height rescaled by the noise dispersion, ctq, af — J \ W{u)\'^u'^'^(Pu 
with a noise filter W{u), = •\/2((Ti/(T2), and 7 = (7i/((To(J2). Equation (2) imphes that the 
number density of the Gaussian noise is a sharply decreasing function of z/, which in turn 
suggests that the noise effect will be almost completely negligible for very high signals u > 5. 

The presence of noise contaminates the SZE sky map, making it obscure to count the 
number density of true signal peaks in many different ways: it changes the amplitudes and 
locations of the signal peaks, and also it creates or compromises signal peaks. To make 
matters worse, the fraction of the signal peaks is usually only a few ten percent of the noise 
peaks. Lee (2002) suggested, however, that all the various noise-contamination effects on 
the number density of real chistcrs nci{y) could be quantified by a one-dimensional Gaussian 
distribution with the unit dispersion p{x) = (l/\/27r)e^^^/^, and that riciiy) be statistically 
estimated by deconvolving p{x) and the zeroth-order approximation to the cluster number 
density n^^\y) = nsz{y) — 'ng{y) where nsz{y) is the number density of local peaks in a SZE 
map: 

Tytot p 

^'S\y)-j;^ p{x)n^{y - x)dx, (3) 

sz J 

where A^*°* and A^*°* are the total number of local peaks in a SZE map and that of real 
clusters, respectively 

The stabilized deconvolution of p{x) and n^^^ (y) can be performed in the Fourier space 

as 

p{yk) 

where yk, Ud, n^^i ^ and p are the Fourier counter parts of n^/, n~^i , and p, respectively. 
Here, Wp{yk) represents the Wiener optimal filter (Press et al. 1992), a key quantity to the 
deconvolution process. It is introduced to stabilize the process of the deconvolution on a 
discretized mesh. The functional form of Wp{yk) is determined from the power spectrum 
of the y-peaks; the ratio of the true power spectrum (P_) to the measured one (-P+). The 
measured power spectrum P+ is usually a corrupted one due to the numerical noise effect. 
The true power spectrum P_ is supposed to be well estimated by subtracting the obvious 
noise tail from the measured power spectrum P+ (see Fig. 2). Finally, the inverse Fourier 
transform of nci{yk) gives the desired estimate of the cluster number density nci{y)- 

The validity and stabihty of the above statistical method depends on the heights of 
y-peaks considering the definitions of n^^\y) and p{x). In the very low y range where the 
noise peaks are dominant (say, y < , where is the noise level), n^^\y) is close to zero, 
and hence the deconvolution likely fails. In other words, in the low y range, most of the local 
peaks can be safely regarded as just noise. Also, in the very high y range (say, y > Scrjv) 



- 5 - 



where the signal peaks are dominant, we have rid {y) ~ ''^sziu) ~ fici{y). In this range the 
noise contamination effect is almost completely negligible so that the dcconvolution of p{x) is 
not meaningful any longer. Hence it is naturally expected that the above statistical method 
works best in the range of medium y peaks (say, (Tn < y ^ ^<^n) where n^^\y) has non-zero 
values while containing a substantial noise-contamination that can be quantified as p{x). 
We emphasize that it is this range of cttv < 2/ < Scxiv where we hope to estimate the cluster 
number density more efficiently than the conventional method. 



3. APPLICATION TO NUMERICAL SIMULATIONS 

We use the SZE cluster catalogues produced from a large cosmological simulation carried 
out by the Virgo consortium (Jenkins et al. 2001; Yoshida et al. 2001). We make maps of 
2048^ pixels in a field of view of 3.14 degree on a side. Details of the generation of the SZE 
cluster maps are found in Yoshida (2002). For our purpose, we use a total of 8 realizations 
of the SZE cluster maps, each of which corresponds to a patch of the SZE sky map. We first 
smooth the SZE cluster maps by a Gaussian filter of scale radius Rg in order to mimic the 
diffuse and extended clusters observed in real surveys. The scale radius Rg corresponds to the 
FWHM (full-width at half maximum) of the clusters which varies with different observations. 
Here we choose Rg = 2 arcminute, given that the natural beam is expected to have a FWHM 
of 2 arcminute in AMiBA experiment. Next, we construct a Gaussian random field on the 
same 2048^ pixels using the Monte-Garlo method, assuming the fiat-sky approximation. To 
simulate the instrumental noise of the real SZE observation, we adopt a white-noise power 
spectrum as the random field, and then convolve it with the noise-cleaning filter Wc which 

— * 

is given in Fourier u-space (u: the Fourier counterpart of9, u= |u|) as (Pen et al. 2002; Lee 
2002): 

1 r / 7,2^2 \ / „2^2_\l 

(5) 



Wciu) = - 
u 



exp I — — exp ' 



where Ob = Q^a with 6a = 0.5 arcminute suitable for AMiBA (Pen 2002, private commu- 
nication). For each of the 8 realizations of SZE cluster maps, wc construct an independent 
noise field, and combine the cluster sources with the noise field to construct a realistic, noise- 
added SZE map. In combining, we rescale each field by the y-dispersion of the noise, . 
Figure 1 shows one realization of the sky maps of original SZE clusters (the top panel) and 
its smoothed version combined with the noise field (the bottom panel), respectively. From 
each realization of the total SZE map, we determine the number density of the peaks Ugz 
as a function of the rescaled peak heights. Using the measured n^^ along with equation (2), 
we compute the zero-th approximation to the cluster number density, n^^\ and measure the 
power spectrum of the peaks from the total SZE map in order to find a Wiener optimal filter 
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Wp. Figure 2 shows an example of the peak power spectrum computed from one reahzation. 
The sohd hne corresponds to the measured power spectrum P+ from the total SZE map re- 
alization, revealing the noise tail. The dashed line represents the estimate of the true power 
spectrum P_, obtained by subtracting the noise tail from P+. The Wiener optimal filter for 
the deconvolution is then approximated as P_/P+. 

Finally, by using this Wiener optimal filter, we perform the stabilized deconvolution of 
the zero-th order signal distribution n^.^"* and the one- dimensional Gaussian distribution p{x) 
to estimate the true number density of the signal peaks. We found that for all realizations 
our method works quite well in the range of ay < y < 5(7y within a ~ 25% error. Figure 
3 plots the number counts of local peaks per degree square as a cumulative function of the 
rescaled Comptonization parameter Nci{> y) = Jynci{y')dy' for the case of each realization 
separately. The long-dashed, the dotted, and the solid lines represent the cumulative number 
densities of local peaks from the total SZE map, noise peaks from equation (2), and the signal 
peaks from the reconstruction, respectively, while the square dots represent the cumulative 
number densities of signals from the original SZE simulations before the noise contamination 
effect added. Figure 4 plots the same as Figure 3 but averaged over all 8 realizations. Our 
method reconstructs also the averaged peak distribution fairly well. 



4. SUMMARY AND CONCLUSIONS 

We have tested the validity of the SZE statistics proposed by Lee (2002) using numerical 
simulations. The statistical method is devised for the interferometric SZE observations in 
drift scanning mode, targeting at AMiBA as a model experiment, and applied to estimate 
the number density of cluster signals that occupy only a small fraction (~ 10%) of the 
noise-contaminated SZE maps more efficiently than the conventional method based on the 
threshold cut-off. The method is developed to model the dominant noise contamination effect 
on the SZE peak distribution by a Gaussian scatter and to deconvolve the noise scatter and 
the SZE peak distribution. The stability of the deconvolution process has been achieved by 
convolution with a Wiener optimal filter estimated from the measurable power spectrum of 
the SZE peaks. 

We have applied the method to a large number of realistic SZE sky maps, and found 
that it indeed works for all realizations, giving a fairly good estimate of the cluster number 
density, including even those clusters whose peak heights are compatible with noise-levels. 
In the conventional approach, those low-amplitude peaks are usually all disregarded as noise, 
so that a large fraction of low-amplitude cluster signals are missing. Our method allows us to 
count those missing low-amplitude clusters as quickly and accurately as possible. Hence, we 
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conclude that it will provide an efficient method to determine the number density of galaxy 
clusters in the future interferometric SZE surveys. 

We thank K. Yoshikawa, U.L. Pen, and Y. Suto for helpful discussions and useful com- 
ments. We both acknowledge gratefully the financial support of the JSPS (Japan Society of 
Promotion of Science) fellowships. This research was funded by the Grant-in-Aid for Scien- 
tific Research of JSPS (02674 and 12640231). The simulation data used in this paper are 
publicly available at www.mpa-garching.mpg.de/galform/virgo/vls 
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SMOOTHING + NOISE 



Fig. 1. — Simulated thermal SZ effect maps. The maps shows the same field of view of 3.14 
degree on a side. Top: a pure thermal SZ effect; Bottom: a smoothed SZE map combined 
with a Gaussian noise field (see §3 in text). 
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Fig. 2. — The power spectrum of the rescaled y-peaks measured from one reahzation of 
the total SZE map. The sohd hne is the measured power spectrum, while the dashed line 
represents our estimate of the true power spectrum. 



- 12 - 



T3 10 

^ 1 

A 1 




0.1 ^ 



100 
^ 10 



A 



1 

0.1 



'o. 100 

10 



A 



0.1 



U 100 

CD 

TS 10 



A 



0.1 



n — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I ^1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



(jN = 4. 87x10-6 
y 

total SZ 

I I I I I I I I I I I 



(7N = 3 73x10-6 
y 



I I f I ii l I I I I I I I I I I I I I 



(jN=3. 72x10-6 

noise 

I I I I I I I I I I I 



I I I I 




(jN=4. 63x10-6 

, signal 

I I I I I I I I I 



I I I I f- I 



(7^=5.54x10-6 
reconst. 



J I I I I I I I I I I I L 



(7^=3.29x10-6 



I I I I I I I I I I I I I 



(7^=3.98x10-6 



I I I I I I I I I I I f- I 



ir (7^=2.99x10-6 



LJ I I I I I I I I I I I I I I I L 



4 



Fig. 3. — The number counts of local peaks per degree square as a cumulative function of 
the rescaled Comptonization parameter in each realization. The long-dashed, the dotted, 
and the solid lines correspond to the number counts of the total peaks of the SZE maps, the 
noise, and the reconstructed signals, respectively, while the square dots correspond to the 
number counts of the true signals from the original SZE map from simulations before the 
combination with the noise field. 
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Fig. 4. — Same as Figure 3 but the averaged over the 8 reahzations. 



